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Abstract 

We report the first observation of two Cabibbo-suppressed decay modes, H+ -^ 
E"'"7r~7r^ and H^ -^ E~7r^7r"'". We observe 59 it 14 over a background of 87, and 
22 lb 8 over a background of 13 events, respectively, for the signals. The data were 
accumulated using the SELEX spectrometer during the 1996-1997 fixed target run 
at Fermilab, chiefly from a 600 GeV/c S~ beam. The branching ratios of the de- 
cays relative to the Cabibbo-favored H+ -^ H^7r+7r+ are measured to be -B(H+ -^ 
E+7r-7r+)/5(H+ -^ H-7r+7r+) = 0.48 ± 0.20, and B{E+ -^ S-7r+7r+)/B(H+ -^ 
H^7r+7r+) = 0.18 it 0.09, respectively. We also report branching ratios for the same 
decay modes of the A^ relative to A^ -^ pK^n^ . 
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1 Introduction 



Studying Cabibbo-suppressed (CS) decays of liadrons provides insights into 
the weak interaction mechanism for non-leptonic decays [1]. Comparing the 
strengths of CS decays to their Cabibbo-favored (CF) analogs, one can, in a 
systematic way, assess the contributions of the various mechanisms. In addi- 
tion, comparing the same or similar decay modes of different baryons allows 
some additional insights. Even though any CS decay mode of the S^ is a CF 
mode of the A^, the detailed arrangement of the different final-state quarks 
into hadrons might be different, as shown in the spectator diagrams in Fig. 1 
and Fig. 2. While in the case of the A+ both final-state baryons, the S"*" and 
the S~, can have the s quark resulting from the CF c decay, in the case of the 
S^ decays with the identical final-state hadrons, only the S~ can be formed 
with the CS c decay product. By comparing several decay modes of different 
hadrons some information about the importances of direct quark emission at 
the decay stage and from quark rearrangement due to final-state scattering 
might be obtained. 
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Fig. 1. Spectator diagrams for the decays (from left to right) H+ — > S^7r+7r+, 
H+ -^ pK~n^, H+ — > E+7r~7r"'", and H+ — > E~7r+7r+. The corresponding 
VF-exchange diagrams and additional final-state quark rearrangements are not 
shown here. 






Fig. 2. Spectator diagrams for the decays A^ — > pK~Tr~^ (left), A+ — > S"'"7r~7r"'" 
(middle), and A^ -^ S~7r"'"7r''' (right). The corresponding 14/^-exchange diagrams 
and additional final-state quark rearrangements are not shown here. 

Modern methods for calculating non-leptonic decay rates of the charm hadrons 
employ heavy quark effective theory and the factorization approximation [2]. 
Nonetheless, the three-body decays of charm baryons are prohibitively diffi- 
cult to calculate due to the complexity of associated final-state interactions. 
Measurements of the relative branching fractions of charm baryon states, both 
CF and CS, give additional information about the structure of the decay am- 
plitude and the validity of the factorization approximation. 

Until now, the only CS S+ decays reported are S^ -^ pK~7i^ [3,4] and S^ — >• 
Tj^K''K~^ [5]. In this letter, we present the first observations of S^ — >■ E^tt^tt^ 
and S^ -^ S^7r"'"7r"'", and determine their branching ratios relative to the CF 
S+ -^ S^7r+7r+. To validate our analysis method, we also report the branching 
ratios 5(A+ -^ ^+tt-tt+)/B{A+ ^ pK-Ti+) andS(A+ -^ ^-tt+tt+)/B{A+ -^ 
S+7r^7r+) and compare them to previously reported results [6-8]. 



2 Experiment 



SELEX is a high energy hadroproduction experiment using a 3-stage spec- 
trometer designed for high acceptance for forward {xp > 0.1) interactions. The 
main goal of the experiment is the study of production and decay properties 
of charm baryons. Particles in the negative (600 GeV/c, ~ 50 % S~, ~ 50 % 



7i~) and positive beam (540 GeV/c, ~ 92 % p, ~ 8 % 7r+) were tagged by a 
beam transition radiation detector. The data were accumulated from a five-foil 
segmented target (2 Cu, 3 C, each separated by 1.5 cm) with a total thickness 
of 5 % of an interaction length for protons. The spectrometer had silicon strip 
detectors to measure the beam and outgoing tracks, giving precision primary 
and secondary vertex reconstruction. Momenta of particles deflected by the 
analyzing magnets were measured by a system of proportional wire chambers 
(PWCs), drift chambers and silicon strip detectors. Momentum resolution for 
a typical 100 GeV/c track was ap/p ~ 0.5 %. Charged particle identification 
was performed with a Ring Imaging Cherenkov detector (RICH) [9], which 
distinguished K"^ from tt"^ up to 165 GeV/c. The proton identification effi- 
ciency was > 95 % above proton threshold (~ 90GeV/c). For pions reaching 
the RICH detector, the total mis-identification probability due to all sources 
of confusion was < 4 %. 

Interactions were selected by a scintillator trigger. The trigger for charm re- 
quired at least 4 charged tracks after the targets as indicated by an interaction 
counter and at least 2 hits in a scintillator hodoscope after the second ana- 
lyzing magnet. It accepted about 1/3 of all inelastic interactions. Triggered 
events were further tested in an on-line computational filter based on down- 
stream tracking and particle identification information. The on-line filter se- 
lected events that had evidence of a secondary vertex from tracks completely 
reconstructed using the forward PWC spectrometer and the vertex silicon. 
This filter reduced the data size by a factor of nearly 8 at a cost of about a 
factor of 2 in charm yield. From a total of 15.2 ■ 10^ interactions during the 
1996-1997 fixed target run about 10^ events were written to tape. A more 
detailed description of the apparatus can be found elsewhere [3,10]. 



3 Data Analysis 



In this analysis, secondary vertex reconstruction was attempted when the x^ 
per degree of freedom for the fit of the ensemble of charged tracks to a sin- 
gle primary vertex exceeded 4. All combinations of tracks were formed for 
secondary vertices (in a first step with x^^^ < 9, but harder cut values were 
applied at later stages) and tested against a reconstruction table that spec- 
ified selection criteria for each charm decay mode. Secondary vertices which 
occurred inside the volume of a target were rejected. Common identification 
criteria for the different decay modes were: proton and kaon candidate tracks 
were required to be identified by the RICH detector to be at least as likely 
as a pion; if a pion candidate track reached the RICH detector, we applied 
as a loose requirement that it had to have a likelihood of at least 10 %, oth- 
erwise it was always accepted; hyperon (S^, S~) decays were identified by 
disappearance of a track in a limited decay interval (5 — 12 m downstream 



from the target), requiring that the candidate track had hits in the tracking 
detectors before the first and in-between the first and second magnet, but no 
hits assigned along the extrapolated trajectory in the 14 chambers after the 
second analyzing magnet; this category of tracks gives unique S'^ identifica- 
tion but is ambiguous between S^ and S^. Additional ambiguities in the mass 
assignments may arise due to loose particle identification criteria for p, K^, 
and TT^; tighter cuts on the identification criteria would reduce the accessible 
momentum range and the number of observed events. 

As additional cuts with variable values depending on the decay modes and 
the relative branching ratio to be determined we used: 

• the separation between the primary and secondary vertices in units of its 
error (L/a) and the error itself (a); 

• the reconstructed charm momentum vector point-back to the primary ver- 
tex, expressed as the square of the distance of the reconstructed charm 
momentum vector to the primary vertex in the target plane in units of its 
error (pvtx); 

• the sum of the squares of the transverse momenta of the daughter tracks 
with respect to the charm hadron direction of flight (Sp^z); 

• the second-largest miss-distance of the daughter tracks in the target plane 
in units of its error (scut); 

• minimum momenta for the tt^ (p^r) and hyperon (phyp) daughter tracks. 

The selection criteria and the actual values for these cuts are discussed in the 
following sections and are listed in tables 2 and 3. 

The total acceptance (geometrical acceptance and reconstruction efficiencies) 
for the different decay modes of interest was estimated by embedding Monte 
Carlo charm decay tracks into data events. Momentum and energy were not 
conserved in the process, but studies indicate this has little effect on the 
single-charm acceptance calculation. Events were generated with an average 
transverse momentum <pt>= l.OGeV/c and longitudinal momentum dis- 
tributions according to (1 — 3;^?)", with n = 2.5 {n = 2.45 ± 0.18 for A+ 
production with a S^ beam [11]). The value of n was varied during the sys- 
tematic studies and did not affect the final branching ratio results. Detector 
hits, including resolution and multiple Coulomb scattering smearing effects, 
produced by these embedded tracks were folded into the hit banks of the un- 
derlying data event. The new ensemble of hits was passed through the SELEX 
off-line software. The acceptance is the ratio of the number of reconstructed 
events to the number of embedded events in a particular mode. For the de- 
termination of the branching ratios only the relative acceptances are relevant, 
leading to a cancellation of most systematic effects associated with the accep- 
tance corrections. 



4 First Observation of Sj" — > S+tt 7r+ and Sj" 
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Fig. 3. Invariant mass distributions of S+tt tt^ (left) and S tt^tt^ (right). 

In fig. 3 we show the invariant mass distributions of S+tt^tt^ and S^7r+7r+, 
over the full mass range evaluated. In each distribution we can see two peaks, 
corresponding to the A+ and S+ decays. The cuts used for the two distributions 
are shown in the first two rows of tables 2 and 3. Additionally we required 
in both channels that at least one of the pions reached the RICH detector. 
For S+7r~7r+ we applied a < 0.10 cm, Xg^^, < 7, phyp > 70 GeV/c and events 
with an invariant mass around the S+ mass in the S+i^~7r+ interpretation 
were removed; for S^7r+7r+ a < 0.08 cm, Xscc < 6, phyp > 80 GeV/c, scut > 
4 and events with an invariant mass around the S+ mass in the S^7r+7r+ 
interpretation were removed. The selection of the individual cut values was 
based on prior SELEX analyzes and tuned to suppress backgrounds mostly 
present in the S^ mass region. To both distributions we adjust the sum of 

Table 1 

Results of the Gaussian parts of the fits to the distributions presented in fig. 3. 



Mode 


Mass [MeV/c2] 


Events 


Mass [MeV/c^] 


Events 


S+7r-7r+ 


2288.1 ±2.2 


74.2 ± 13.8 


2471.6 ±3.9 


58.7 ±13.5 


S-7r+7r+ 


2286.0 ± 1.8 


46.4 ±10.1 


2463.3 ± 3.0 


22.3 ±7.5 



two Gaussians with fixed widths (given by Monte Carlo) and a second degree 
polynomial. The results of the fits are summarized in table 1. For the decay 
S+ -^ S+7r^7r+ we observe S = 58.7 ± 13.5 signal events over a background of 
B = 87.3±6.7, corresponding to a significance S/^/B = 6.3±1.5. For the decay 
S+ -^ S^7r+7r+ we observe S = 22.3 ± 7.5 signal events over a background 
oi B = 12.8 ± 2.5, corresponding to a significance S/yB = 6.2 ± 2.2. The 
masses of both the A+ and the S+ are slightly higher (in the case of S+tt^tt^) 



and lower (E~7r"'"7r"'") than the nominal values. Varying the bin width and the 
fixed widths of the Gaussians also gives consistent results. We verified that the 
number of observed events varies as a function of the cut variables, especially 
L and L/a^ in the same way as expected from Monte Carlo (e.g. the observed 
events have the lifetime of the Sj"). 



5 Measurement of Branching Ratios 



For the different branching ratio measurements we used different selection cuts, 
chosen under the criteria to minimize systematic effects on the final result. We 
selected central cut values within a region where the Monte Carlo described 
well the distributions of all the observables, for both the decay mode of interest 
and the normalization mode. If the same mode was used in different branching 
ratio determinations, this selection could result in a different set of cuts. We 
only used interactions initiated by a S^ as beam particle and all decay prod- 
ucts (with the exception of the hyperons) had to be within the RICH accep- 
tance and correspondingly identified. In the different modes we removed events 
stemming from the following reflections due to ambiguities in the mass assign- 
ments: D+ -^ K~n+Ti+ (1), D+ -^ K+K-Tr+ (2), D+ ^ K+Tr~TT+ (3), Df -^ 
K+K-TI+ (4), Dt -^ K+71-71+ (5), A+ -^ p^-^+ (6), A+ -^ S-7r+7r+ (7), 
A+ -^ pK'n+ (8), S+ -^ S"7r+7r+ (9), S+ -^ ^+K-ti+ (10). We indicate the 
removed reflections and the corrections due to the removal in tables 2 and 3. 

In flg. 4 we show the invariant mass distributions of S+tt^tt^, S^tt+tt^, 
pK^Tf^, and S^7r+7r+ in the S+ mass region. 

To verify our analysis method we determine the relative branching ratios of 
two A^ decay modes which are identical to our newly observed S^ modes, 
using the fact that every CF A+ decay mode is also a CS S^ mode. In flg. 5 
we show the invariant mass distributions of S+7r^7r+, S^7r+7r+, and pK^n^ 
in the A^ mass region. 

The number of observed events was determined by adjusting a Gaussian of 
flxed width (given by Monte Carlo) and a flrst-order polynomial to the dis- 
tributions shown in flgs. 4 and 5 ^^ . While most of the reflections lie outside 
of the mass peaks and are removed to smoothen the backgrounds, some of 
them extend below the peaks and remove good events; we studied this effect 
carefully with Monte Carlo and corrected the number of observed events for 
these losses. To determine the correction we simulated the shape of the invari- 



^^ Counting the number of entries above the extrapolated background, and using 
a second-order polynomial for the background, we obtain within errors the same 
number of events. 
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Fig. 4. Eight invariant mass distributions of: pK^n^, S+7r^7r+, S^7r+7r+, 
H^TT^Tr^, used to determine the four relative branching ratios (in 
pairs from top to bottom) B{E,'^ -^ pX^7r+)/i?(S+ -^ H^7r+7r+), 
B{E+ -^ S+7r-7r+)/5(H+ -^ E;-7r+7r+), B{E+ -^ S-7r+7r+)/5(H+ -^ H-7r+7r+), 
and B{Ef -^ S^7r^"7r+)/-B(H+ -^ S+7r~7r"'"), respectively. Different selection cuts 
were used for each branching ratio (see text). We adjust a Gaussian (fixed width 
given by Monte Carlo) over a linear background to each of the distributions. The 
event yields are summarized in table 2. 



ant mass distribution of the reflected mode (including all cuts) and scaled the 
number of events below the peak region to the number of observed events in 
the reflected mode, keeping the same relative error for the number of observed 
events. We also studied correlations for the cases where more than one reflec- 
tion was removed and found them to be small and negligible. The corrected 
yields for the different modes are, along with the cuts used to obtain the distri- 
butions and the corresponding total acceptances, presented in tables 2 and 3. 
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Fig. 5. Six invariant mass distributions of: S+vr 7r+, pK 7r+, S 7r+7r+, used 
to determine the three relative branching ratios (in pairs from top to bottom) 



B{Ai 



S+ 



vr vr 



-)/BiAi 



pK- 



vr ' 



B{A+ -^ T,-TT+TT+)/B{Aj -^ S+7r-7r+), 
and B{A^ -^ T,~tt'^tt~^)/B{A'^ -^ pK~iT~^), respectively. Different selection cuts 
were used for each branching ratio (see text). We adjust a Gaussian (fixed width 
given by Monte Carlo) over a linear background to each of the distributions. The 
event yields are summarized in table 3. 

To obtain the branching ratios, we divided the number of observed (corrected) 
events of the two modes, and divided again by the relative acceptance. The 
statistical error on the acceptance is negligible, and most systematic errors 
cancel in the relative acceptance. 

For the systematic studies we varied any single cut value, as well as the pa- 
rameter n for the xp distribution in the Monte Carlo simulation, within some 
range and determined the branching ratio for every set of parameters; for the 
set of cuts used we did not observe evidence of any trend; all systematic vari- 
ations are small compared to the statistical error and will be ignored in the 
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Table 2 

Number of observed events and total acceptances for the different H+ decay modes, 
with the corresponding cuts applied to each mode. Common cuts are: Xsec < 8' 
a < 0.10 cm, scut > 8, Phyp > 40GeV/c. The first two rows refer to the signals 
shown in fig. 3, with different common cuts as described in section 4. "Corrected 
Events" are the number of observed events plus the corrections due to the removal 
of reflections and are shown separately in parenthesis; we keep the relative error 
from the fits. 



Mode 


L/a 


pvtx 


[GeVVc2] 


Removed 
Reflections 


Corrected 
Events 


Acceptance 

[%] 


^ + TT-TT+ 


> 12 


< 13 


>0.4 


(10) 


58.7 ±13.5 


0.750 


s-^+^+ 


>8 


< 10 


>0.5 


(9) 


22.3 ±7.5 


0.950 


pK^TT^ 


> 11 


< 13 


>0.3 


(1,2,4,6) 
(7) 


(47.4 ± 14.0) ± 14.0 
(67.0 + 2.4) ± 10.9 


4.164 
0.915 


S+7r-7r+ 
S-7r+7r+ 


> 13 


< 13 


>0.35 


(6) 

(7) 


(20.7 + 1.6) ±8.6 
(63.3 + 2.0) ±10.4 


0.586 
0.825 


H-7r+7r+ 


> 13 


< 10 


>0.35 


(2,3,5,9) 

(7) 


(9.5 + 5.0) ±6.4 
(61.1 + 2.6) ±9.3 


0.988 
0.800 




> 13 


< 10 


>0.35 


(2,3,5,9) 
(6) 


(9.5 + 5.0) ±6.4 

(18.4 + 1.5) ±7.5 


0.988 
0.570 



final results since they would not affect the quadrature sum of the total error. 

The resulting branching ratios are shown, together with previously measured 
values, in table 4. 
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Table 3 

Number of observed events and total acceptances for the different A+ decay modes, 
with the corresponding cuts applied to each mode. Common cuts are: Xsec < 4' 
G < 0.10 cm, scut > 8, phyp > 40GeV/c. The first two rows refer to the signals 
shown in fig. 3, with different common cuts as described in section 4. "Corrected 
Events" are the number of observed events plus the corrections due to the removal 
of reflections and are shown separately in parenthesis; we keep the relative error 
from the fits. 



A+ 
Mode 


L/a 


pvtx 


Spi2 

[GeVVc^] 


Removed 
Reflections 


Corrected 
Events 


Acceptance 

[%] 


E+7r-^+ 


> 12 


< 13 


>0.4 


(10) 


74.2 ± 13.8 


0.450 


s-^+^+ 


>8 


< 10 


>0.5 


(9) 


46.4 ±10.1 


0.500 


pK^TT^ 


> 11 


<7 


>0.3 


(8) 


(46.6 + 3.0) ± 9.6 
(561.2 + 0.0) ±25.5 


0.292 
2.367 


pK-7r+ 


> 11 


<4 


>0.4 


(5,9) 


(29.7 + 2.2) ±6.6 
(450.7 + 0.0) ±22.3 


0.434 
1.923 


S-7r+7r+ 


> 11 


<4 


>0.4 


(5,9) 
(8) 


(29.7 + 2.2) ±6.6 
(43.0 + 3.4) ±8.1 


0.434 
0.241 
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6 Discussion and Conclusions 



In table 4 we summarize the results for the different branching ratios measured 
in this work. Comparing our results with previously measured ones (where 
available) shows good agreement. 

Table 4 

Results of the different Branching Ratios measured in this analysis, and comparison 
to previously published results (if available). Also shown is the a-parameter (see 
text) for each branching ratio result. 



Branching Ratio 


This Analysis 


Other Measurements 


S(S+- 


S+7r-7r+) / 


0.48 ± 0.20 


- 


BiEt - 


■> H^7r+7r+) 


a = 6.4 ±2.7 




BiEt^ 


S-7r+7r+) / 


0.18 ±0.09 


- 


B{Et - 


^H-7r+7r+) 


a = 2.5 ±1.2 




B{Et^ 


S-7r+7r+) / 


0.42 ± 0.24 


- 


B{Et- 


^ S+7r-7r+) 


a = 0.43 ±0.25 




B{Et- 


>pK-TT+) 1 


0.194 ±0.054 


0.234 ± 0.047 ± 0.022 [4] 


B{Et 


-^ H-7r+7r+) 


a = 2.6 ±0.7 


0.20 ± 0.04 ± 0.02 [3] 


B{K+ ^ 


s-^+^+) / 


0.314 ±0.067 


- 


B{Kt 


^pK-^+) 


a = 0.30 ±0.07 




B{kt - 


s+^-^+) / 


0.72 ±0.14 


0.74 ± 0.07 ± 0.09 [7] 


S(A+ 


^pK-TT+) 


a = 0.68 ±0.14 


OMt'oll [6] 


S(A+ ^ 


s-^+^+) / 


0.38 ±0.10 


0.53 ±0.15 ±0.07 [8] 


B(A+ 


^ S+7r-7r+) 


a = 0.39 ±0.11 





To quantify the effects of final-state quark rearrangements in the different 
decays via the relevant relative matrix elements, we calculate a, which is 
defined as the measured relative branching ratio corrected for phase space 
differences and, in the case of comparing CF and CS modes, for the ratio of 
the CKM matrix elements (Kd/Ks = 0.233 ± 0.001 [12]). We note that the 
a-parameter for i?(E~7r"'"7r"'")/-B(S"'"7r~7r"'") is consistent in the decays of both 
the A^ and the S+. Comparing the diagrams in fig. 1 and 2 we conclude that 
the source of the final-state quark does not affect the relative matrix element 
significantly. 



In summary, we observe for the first time the Cabibbo-suppressed decay modes 

S^TT+TT^ and estimate their branching ratios. With 



S+ -^ S+TT 7r+ and S^ 
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the same analysis method we also analyze previously reported modes of both 
the S^ and the A+ and find good agreement. 
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